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> ; Abstract 



The 't Hooft model (two-dimensional QCD in the limit of large number of 



<N 
in 

<^> . colors) is used as a laboratory for exploring various aspects of the heavy quark 

expansions in the nonleptonic and semileptonic decays of heavy flavors. We 
perform a complete operator analysis and construct the operator product ex- 
Q_; pansion (OPE) up to terms 0(l/mg), inclusively. The OPE-based predictions 

for the inclusive widths are then confronted with the "phenomenological" re- 
sults, obtained by summation of all open exclusive decay channels, one by 
one. The summation is carried out analytically, by virtue of the 't Hooft 
equation. The two alternative expressions for the total widths match. We 
^ ' comment on the recent claim in the literature of a 1/mg correction to the 

total width which would be in clear conflict with the OPE result. 

The issue of duality violations both in the simplified setting of the 't Hooft 
model and in actual QCD is discussed. The amplitude of oscillating terms is 
estimated. 
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I. OVERVIEW 



The development of heavy quark theory started in the 1980's, has essentially been com- 
pleted. While at the early stages the main emphasis was placed on the symmetry aspects 
(the so-called heavy quark symmetry), the present (mature) stage deals with dynamical 
aspects. A formalism based on Wilson's operator product expansion (OPE) [|IJ has been 
developed and applied to many cases of practical interest, in particular to inclusive decays 
of heavy flavor hadrons. The theory of such decays is at a rather advanced stage now (see 
and references therein). Calculations we could not even dream of several years ago have 
become possible. 

The decays of heavy flavor hadrons Hq are shaped by nonperturbative dynamics. While 
QCD at large distances is not yet solved, considerable progress has been achieved in this 
problem. The width of an inclusive transition Hq — > f is expressed through an OPE. 
The nonperturbative effects are then parameterized through expectation values of various 
local operators Oi built from the quark and/or gluon fields. Observable quantities, such as 
semileptonic and nonleptonic widths of heavy hadrons Hq, are then given by 

where Cj are the OPE coefficients, and /i stands for a normalization point separating out soft 
contributions (which are lumped into the matrix elements (ifg|(9j(//)|f7g)) from the hard 
ones (which belong to the coefficient functions q). 

There are many subtle and interrelated issues, both conceptual and technical, associated 
with the operator product expansion in QCD: 

1. Eq. (|1|) represents an expansion in powers of 1/mg with tuq being the Q quark mass 
and the coefficients Cj scaling like (I/wq)*" 2 for an operator with dimension <ij. In 
Th q there are two sources for contributions depending on powers of 1/mg, namely 
higher- dimensional operators and higher-order terms in the expansion of their ex- 
pectation values. In addition every coefficient is a series in the running coupling 
a s( m Q) l/l°g m g) of which only a few terms are known for a given coefficient Cj. 
This immediately raises a grave concern: how can we retain terms suppressed by pow- 
ers of 1/mg without a complete summation of the parametrically larger powers of 
1/ log tuq in the leading coefficient? 

2. Although the normalization point /i conceptually represents a straightforward "book- 
keeping" device for separating hard and soft contributions, it is technically difficult 
to actually carry out such a program since no user-friendly definition of what is soft 
and hard exists in QCD. So far, the vast majority of all discussions related to the 
introduction of /i are conducted in a hand-waving manner. 

3. It is quite conceivable that there are hard nonperturbative contributions in the co- 
efficient functions: c ™ onpert ^ (AqcDA^g) 15 with S being some positive number. The 
possible size of such contributions is essentially unknown. A related problem is the 
convergence (or divergence) of the perturbative series for the coefficient functions Cj. 
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4. Truncating the series (|T|) at some finite order introduces an error estimated by so- 
called exponential terms, which in Euclidean domain look as expressions of the type 
exp[— (tbq/Aqcd)* 1 ]. In order to obtain Y Hq we analytically continue from the Eu- 
clidean domain, where the OPE is well defined and the coefficients q are real, to the 
Minkowski domain where they acquire an imaginary part. Such analytic continuation 
is implicit in Eq. ([I]) and is based on the assumption of smoothness. Under analytic 
continuation the exponential terms convert themselves into oscillating terms of the 
type cos[(mQ/Aq CD ) fc ] 0; the expansion ([|) does not account for them. It can thus be 
understood on general grounds that duality violation is described - or at least modeled 
- by oscillating expressions. To which degree those are suppressed by powers of 1/ rriQ 
depends on details of the strong interactions and the specifics of the process. 

All these questions are circumvented in the so-called practical version of OPE [|J rou- 
tinely used so far in all instances when there is need in numerical predictions. This version 
is admittedly approximate, however. The questions formulated above are legitimate; they 
deserve to attract theorists' attention, and continue to cause confusion in the literature. 
They have to be addressed also because they are emerging as a major source of the un- 
certainties in quantitative predictions; these problems have specifically been suspected to 
underlie phenomenological difficulties encountered recently, e.g. a relatively short lifetime 
of beauty baryons and a relatively small semileptonic branching ratio of beauty mesons. 

We find it useful and instructive to study all these issues in models that while retaining 
basic features of QCD - most notably quark confinement - are simpler without being trivial 
and can be solved dynamically. QCD defined in one time and one space dimension - hereafter 
referred to as 1+1 QCD - is especially suitable for this purpose: with the Coulomb potential 
necessarily growing linearly in two dimensions, quark confinement is built in. Likewise the 
theory is superrenormalizable, i.e. very simple in the ultraviolet domain. There are no 
logarithmically divergent "tails" in the Feynman graphs. As a result, the book-keeping of 
OPE (separation of the hard and soft parts) becomes simple, and all subtle aspects in the 
construction of the OPE can be studied in a transparent environment. 

In particular, the perturbative contributions in the coefficients Cj become an expansion 
in g 2 /rn 2 Q (where g is the gauge coupling in 1+1 QCD). They are thus power-suppressed 
in the same way as the higher-dimensional operators; the first problem formulated above 
therefore does not arise here. Without the logarithmic UV tails the second problem becomes 
tractable. Concerning the third problem it is easy to see that in 1+1 QCD nonperturbative 
corrections cannot generate power suppressed terms in the coefficients q. For the leading 
operator QQ we will find its coefficient function to all orders of perturbation theory (in 
the limit of N c — > oo), demonstrating the convergence of the perturbative series. At the 
same time, the divergence of the condensate expansion in high orders will become manifest 
indirectly, through the occurrence of oscillating terms in Th q , which appear with suppression 
factor (1/mg) 7 in the case at hand. Thus all the four problems formulated above will be 
answered! 

We will perform our explicit calculations for 1+1 QCD in the limit of a large number of 
colors N c - the famous 't Hooft model 0-0. For N c — > oo only planar diagrams contribute 
in QCD; 1+1 QCD has the additional special feature that one can choose a gauge such that 
there are no gluon self-interactions. Then only planar ladder diagrams have to be considered, 
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and we have an exactly solvable theory in our hands. All hadronic matrix elements of 
interest are therefore calculable. This enables us to describe every given transition in two 
complementary ways: we can confront the OPE-based expression with a "phenomenological" 
representation for the same process obtained by saturating the rate by exclusive hadronic 
channels. 

We want to take advantage of these unique features of the 't Hooft model to illustrate 
all crucial elements of heavy quark theory and the theory of inclusive heavy flavor decays in 
particular. One should keep in mind that heavy quark theory, as we know it now, is merely 
an adaptation of the general OPE-based approach. Some of the questions to be discussed 
below can therefore be actually formulated in a wider setting. 

The 't Hooft model has been exploited as a theoretical laboratory for testing various 
analytic QCD methods in applied problems before. Heavy quark symmetry and heavy 
flavor decays were analyzed in Refs. jSHlQfl- The model was used recently for discussing 
general aspects of OPE (convergence of the OPE series, exponential terms violating duality, 
and so on) |TT| , |12| • 



In Ref. [^OJ heavy flavor inclusive widths were calculated numerically, by adding the ex- 
clusive channels one by one. It was found that the inclusive width Th q approaches its asymp- 
totic (partonic) value, and the sum over the exclusive hadronic states converges rapidly. At 
the same time, small deviations from the asymptotic value observed in the numerical analy- 
sis [|l(J were claimed to be a signal of 1/ itlq corrections in the total width, in contradiction 
with the OPE-based result. 

In this work we treat the very same problem, inclusive heavy flavor decays in 1+1 QCD, 
analytically. We first develop a technique perfectly parallel to that in four-dimensional QCD 
0. It includes such elements as a complete operator analysis and the construction of the 
transition operator. Unlike four- dimensional QCD, the coefficient functions for the leading 
operator are exactly calculable (in the limit N c — > oo). Moreover, all relevant expectation 
values of the local operators involved in the problem are calculable too. We get a complete 
prediction through order I/tjiq. 

Then we carry out a "hadronic calculation" of the same width, by saturating all open 
decay modes, using the 't Hooft equation [Q. By comparing the phenomenological represen- 
tation of the total width with the OPE-based formula, we are able to identify, term-by-term, 
the subsequent terms of the heavy quark expansion. The situation actually turns out to be 
simpler than one could expect a priori: 

• In the 1/rriQ expansion for the inclusive width corrections of the order (1/mg) 2 , 
(1/mg) 3 and (I/ttiq) 4 to the parton width come only from the leading operator 
QQ, i.e., from the expansions of its OPE coefficient cqq and its expectation value 
(Hq\QQ\Hq) . Operators of higher dimension contribute to the total width first at 
order (1/mg) 5 . 

• The perturbative series in g 2 /trig for the OPE coefficient of the operator QQ is com- 
pletely defined by the one-loop renormalization of heavy quark mass. The result can 
be formulated in terms of the light-cone gauge formalism as the absence of renormal- 
ization. 

These results are based on a general operator analysis. On the "phenomenological" side, we 
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use a sum rule, which is a consequence of the 't Hooft equation, to show that the total width 
is determined by a quantity coinciding with the matrix element of the cqqQQ term in the 
OPE expansion ([I]), through order (1/mg) 4 . Thus, we observe a perfect match between the 
expression derived from the OPE and from adding up all relevant hadronic channels up to 
high order power corrections. 

After testing the validity of OPE, we exploit results obtained en route in order to discuss 
the issue of oscillating contributions related to the high-order tails in the OPE series that 
are factorially divergent. Due to the simplicity of the model we can estimate them reliably. 
A non-monotonous duality-violating component of the width for large tuq is suppressed by 
high power of 1/toq which we determined. Implications of our analysis for real QCD are 
briefly discussed. 

The remainder of the paper is organized as follows: after formulating the problem in 
Sec. [n] we construct the OPE and calculate the coefficients in Sec. [II; after establishing the 
match between the OPE-based result for the inclusive width and the sum rules for the same 
width resulting from the 't Hooft equation through order (1/mg) 4 in Sec. [TV], we discuss an 
appearance of oscillating terms in the order (1/mg) 7 and the duality violations they cause 
in Sec. [V]; in the same section we discussed along similar lines a possible pattern of the 
violation of the local duality for r decays in 1+3 dimensions; in Sec. [VJ we comment on the 
paper |I| and analyze effects due to nonvanishing masses of light quarks; Sec. |V11| presents 
a general discussion and conclusions. 



II. PRELIMINARIES 

We start by formulating the problem and introducing our notation and conventions. 
In two-dimensional QCD the Lagrangian looks superficially the same as in four dimen- 
sions 

A+i = ~ G%G% + TP - mttyi , iD, = id, + A a ^T a ; (2) 

T a denote generators of SU(N C ) in the fundamental representation, the gluon field 
strength tensor and ipi the quark field (i is a flavor index) with a mass ra^; g the gauge 
coupling constant. 

One has to keep the following peculiarities in mind: g carries dimension of mass as does 
ifjifj. The field strength G a on the other hand has dimension M 2 in our normalization, 
just as in four-dimensional QCD. With the theory being superrenormalizable no (infinite) 
renormalization is needed; observables like the total width Th q can be expressed in terms 
of the bare masses mi and bare coupling g appearing in the Lagrangian. Anticipating the 
large N c limit we will use a parameter (3 instead of g where 

j3 2 = — ( N c - —) with lim (3 2 = finite . (3) 

This dimensionful quantity /3, which - in contrast to m, - provides an intrinsic mass unit 
for the 't Hooft model, can be seen as the analog of Aqcd of four- dimensional QCD. 
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We need at least two quarks denoted by Q and q with masses rriQ and m q , respectively, 
to realize heavy flavor transitions Q — > q. For quark masses we impose 

m Q - m q > , (4) 

where both m q ^ and m q = are allowed for. Condition (f|) guarantees that the inclusive 
methods of Ref. p[ are applicable since it makes the energy release [] in the weak decay 
large relative to the intrinsic scale (3. We will actually employ the dimensionless ratio P/tuq 
as our expansion parameter. 

Next we need to introduce a flavor-changing weak interaction; we choose it to be of the 
current-current form: 

Cak = - J= {TtyQ) fa-fti) ■ (5) 

Here G is an analog of the Fermi coupling constant; it is dimensionless in two dimensions. 
The fields tp at b can be either the light quark or the lepton fields to describe nonleptonic or 
semileptonic decays, respectively. In 1+1 dimensions the axial current reduces to the vector 
one. The most general current-current interaction contains an additional term where the 
vector currents are contracted via the antisymmetric instead of g^ u . For the total width 
- our main focus here - such an additional term is of no importance. The product of scalar 
densities, on the other hand, is inequivalent to that of vector densities; we will briefly discuss 
it, but mainly focus on the V x V interaction (|5|). 

For N c — > oo factorization holds; i.e., the transition amplitude can be written as the 
product of matrix elements of the currents qj^Q and Vyy^Vv F° r the inclusive widths 
which are discussed below the property of factorization can be expressed as follows: 

M Hq T Hq =lmj d 2 xt (H Q \ T {£ wcak (x)4 cak (0)} \H Q ) = G 2 J d 2 xImIV(x) ImT^(x) (6) 

where H^x) and T^ix) are defined as 

rVO) = i (0|T {$ a {x)j^ b {x) ^(0)7^(0)} |0> , (7) 

T^(x) = i{H Q \T{q(x)rQ(x)Q(0) lv q(0)}\H Q ) . (8) 

This factorization follows from the fact that at N c — > oo there is no communication between 
i'al^b and qj^Q currents: any gluon exchange brings in a suppression factor l/N 2 . 

The only difference between the semileptonic and nonleptonic widths resides in Yl^^x), 
in the first case ip a are leptonic fields while in the second case they are the quark fields. 
At = we get one and the same (x) (up to the overall normalization factor iV c ) as 
we will show shortly. For this reason at = the distinction between the nonleptonic 



It can hardly be overemphasized that it is the size of the energy release rather than of tuq that 
controls the reliability of the expansion in four dimensions as well. 
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and semileptonic cases is actually immaterial. At ^ quark and lepton polarization 
tensors U^ u (x) become different. This difference is proportional to powers of m^. It will be 
discussed in Sec. |VT[ For the time being we will treat ip's as massless leptons. 
The one-loop graph determining IL^ is depicted in Fig. [I]. 




FIG. 1. Polarization operator for lepton current 



For a massless fermion ip we get the well-known expression, 

IV(g) = J d 2 xe^IV(x) = ~ - . (9) 

This expression obtained from a one- loop graph is known to be exact. If ip is a lepton 
field, this statement is trivial. If ip is the quark field all gluon insertions inside the loop 
automatically vanish due to special properties of the two-dimensional 7 matrices 0. Thus, at 

= the only distinction between ipa^ being quark rather than lepton fields is an overall 
factor N c on the right-hand side of Eq. ([]) . 

A remarkable feature of Eq. (^) is the occurrence of the pole at q 2 = 0, which is specific 
for the vector interaction. This means that a pair of massless leptons produced by the vector 
current is equivalent to one massless boson, whose coupling is proportional to its momentum 
q^. In the case of the quark fields, it is known || from the early days of the 't Hooft model 
that the vector current ipal^ipb produces from the vacuum only one massless meson, the 
pion. This is readily seen by inspecting the 't Hooft equation ||. 

For all computational purposes the vector current ^al^^b in Eq. (||) can thus be substi- 
tuted by d u (j) / y/n where denotes a pseudoscalar massless noninteracting field, 

^eak = -^^Q^0. (10) 

In other words, the problem is formulated as the inclusive decay of the heavy quark Q into 
a lighter quark q plus a sterile boson (p. More exactly, we deal with the decays of a Q 
containing hadron Hq into a q containing final hadronic state X q plus </>. 

Let us pause here for two remarks, (i) The fact that the interaction vertex of the 
massless field (f) involves e^, see Eq. (0), is most obvious when ip a ^ are quark fields. For in 



2 Namely, one uses the fact that 7 a 7 M 7 a = and any odd number of 7 matrices reduces to one. 
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this case is the pion, as mentioned above, and the pion is coupled to the vector current 
i'aj^b obviously through & v . The case of the leptonic fields ipa,b is indistinguishable; 
therefore, the coupling of is the same, (ii) To keep the analysis to be presented below as 
clean and transparent as possible we want to be free of annihilation and Pauli interference 
contributions to the total width (at least through 0(l/mg)). This is readily achieved by 
assuming throughout the paper that the spectator light quark q sp in Hq is distinct from q. 

In the leading approximation the transition operator is determined by the diagram of 
Fig. where the wavy line corresponds to the <fi quantum. A straightforward calculation 
yields for the transition operator 



T = c° 5o QQ; 2Imc° 



r 



Q 



Aix 



m 



Q 



m Q 



(11) 



where Tq is the decay width for a free quark Q as evaluated in the parton model. 



Q 




Q 



FIG. 2. Transition operator in the leading order. The wavy line of massless <fi field substitutes 
the propagation of lepton pair 



This parton expression will serve as reference in analyzing the (l/mo) n corrections to 
the total width Y. 



III. OPERATOR PRODUCT EXPANSION FOR INCLUSIVE WIDTHS 

A. Catalogue of operators 

The 1 / rriQ expansion for inclusive widths of heavy flavor hadrons is constructed from the 
Lorentz invariant weak transition operator ]13| 



f(Q -Q) = J d 2 x iT {£ wcak (a;)4 eak (0)} = £ a^Ofo) ■ (12) 

The local operators Oi are ordered according to their dimensions. The leading one is QQ 
with dimension cLqq = 1. Higher operators have dimensions di > 1. By dimensional counting 
the corresponding coefficients are proportional to (l/mg)^ -2 ). The ratio of the coefficients 

Ci/cQq is proportional to mQ l+dQQ . 

The coefficients Cj are determined in perturbation theory as a series in f3 2 /m^. It 
is crucial that these coefficients are saturated by the domain of virtual momenta ~ niq 
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and are infrared stable by construction. (All infrared contributions reside in the ma- 
trix elements of the operators 0{.) At this point we should mention a drastic distinc- 
tion between four- and two-dimensional QCD. In four dimensions the expansion parame- 
ter for the coefficients is the running coupling a s (mg); nonperturbative contributions to 
the coefficients coming from distances ~ l/ m Q could in principle show up in the form 
exp(— C / a s {mQ)) ~ (Aqcd/^q) 5 where 5 is some unknown positive index, not necessarily 
integer. In two-dimensional QCD such terms cannot appear: an analog of the exponential 
term above would be exp(— Cmq/ f3 2 ). 

A note concerning the choice of the normalization point \i: we will imply that 

m Q > n > . (13) 

In this range there is no real dependence on /i, so we will suppress the argument \x both in 
the coefficient functions and operators. 

The coefficient functions are not the only source of a itlq dependence. The matrix 
elements of the operators Oi contain an implicit tuq dependence too. (We recall that in our 
formalism, unlike HQET |14 1 , the fields of which the operators Oi are built are the standard 



Heisenberg operators, rather than asymptotic in tuq.) In particular, for the leading operator 
QQ we have the following relation: 



d 2 x QQ = d 2 x { Q 7 oQ + Q Y\ v Q , (14) 




where 7r M = iD^ — niq and the integration over x allows us to omit terms which are total 
derivatives. 

In the rest frame of the hadron Hq the expectation value of QjoQ counts the number of 
Q quarks, 

-L-{H Q \Q>toQ\H Q ) = l. (15) 
zm h q 

The factor 1/2Mh q will be present in all matrix elements; it corresponds to a relativistic 
normalization of the states, 

(H Q (f)\H Q (p)) = 2E HQ 5(f - p) . (16) 

From relation (0) the matrix element of QQ is therefore unity, up to a quadratic cor- 
rection: 

' 'H Q \QQ\H Q ) = l + o(\\ . (17) 



Moreover relation fllil) provides an operator form for 1/ttIq corrections. They come from 
Qit^ir^Q /2iriQ which equals QD\Q/2mn up to I/ttlq corrections. Notice, that the operator 
QD\Q is Lorentz noncovariant and cannot enter directly into the OPE for the total width 
but, as we see, enters indirectly through the matrix element of the operator QQ. 
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Eq. (|14]) contains also the 'chromomagnetic' operator and it looks as this operator con- 
tributes to I/tuq corrections. It is not, however, the case. Indeed, this operator can be 
rewritten as follows 

O g = ~ Qcr^G^Q = % - Q^G^Q = % - Q^G^l - 7 °)Q + ~ 0(1 - i'H^G^Q , 

(18) 

where the relation a^ v = ^e^ v is used. Taking advantage of the non-relativistic equations 
of motion to replace 

1 



j°)Q by (l/mQ)j 1 iDiQ we get (up to total derivatives) 



Or 



Irrir 



Q(D»G^)YQ + O 



m, 



2m, 



m, 



(19) 



where t a stand for the generators of the color group SU(N C ). Thus, the operator Oq reduces 
to a four-fermion operator 0± q with coefficient g 2 /m,Q. 

The absence of operators with the gluon field strength tensor G^ v in the OPE is a 
specific feature of two-dimensional QCD. The physical reason for the reducibility of the 
gluonic operators is the absence of real gluons in two dimensions. A particular consequence 
of Eq. (|i9|) is that in Eq. ( |H]) the chromomagnetic operator generates 1 / m?q terms only. 

Thus we come to the following representation for the matrix element of the leading 
operator QQ, 



[H Q \QQ\H Q ) 
2M Ho 



1 (Hq\Q(-DI)Q\H q ) 



2mj 



2M Hc 



+ 



g 2 (H Q \Q lt ,t a QE q q rt a q\H Q ) | / 1 



2vP? r 



2M, 



(20) 



Let us proceed further with the operator analysis. The first subleading operator is the 
dimension-two four-fermion operator of the type 



Aq = QT 1 QqT 2 q 



(21) 



where r 12 denote color and spinor matrices. This is in distinction with 1+3 QCD where 
the first subleading operator was Qa^G^Q. On dimensional grounds the operator (9 4g , if 
present, could produce a linear 1/mq correction in the total decay width. To this end the 
corresponding coefficient must arise in zeroth order in the coupling g 2 (see the diagram in 
Fig. I). 




FIG. 3. Four-fermion operators in the leading order 
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For the V x V weak coupling of Lagrangian (|T0D this graph vanishes identically in the 
imaginary part, Im c iq = 0, provided the leptons are massless, see discussion in Sec. ||. (If 
the fields propagating in the loop are massive, the operator 0± q appears. The corresponding 
modifications are considered in Sec. [V|.) 

However, in the case of a scalar-scalar weak interaction of the form 

^weak = (qQ) Mb) (22) 

the graph of Fig. is not zero, and gives the following contribution to the transition operator: 

ImT s = C f(Qq)(qQ). (23) 

With this contribution the total width takes the form (we put m q = for simplicity): 

r s GlmQ L (H Q \(Qq)(qQ)\H Q } \ 
Hq 16vr ^ + m Q (Hq\QQ\H q ) ) ■ { ] 

Let us emphasize that, unlike four- dimensional QCD where no operator can induce a 



1/rn.q contribution to the total width |15[,16l, this can happen in two dimensions. The 



vanishing of Im C4 q for the V x V weak Lagrangian in the leading order is a specific dynamical 
feature of this particular Lorentz structure of the weak interaction. Note at this point that 
the argumentation presented in Ref. |l5j was not sufficient to prove the absence of 1/mg 



corrections in real QCD, see the discussion in Sec. |V A| . 

In Sec. [T UB 3| we will show that Im c^ q vanishes not only in the leading order in strong 
coupling but also in the order g 2 . The first nonvanishing contribution to Im c^ q comes in 
the order g A what leads to I/tuq corrections to the width. 

Next in the list comes the dimension-three operator containing six quark fields: 

{QT 1 Q)(q a T 2 qb){qc^^qd) ■ 

Along the same line of reasoning as for four-fermion operators we show that six-fermion 
operators appear only in the order g 6 . It therefore contributes in order l/m%. For multi- 
fermion operators every extra qq pair brings in an extra 1/mg suppression. 

In summary: to a quite high accuracy the operator QQ is the only one to contribute. 



B. Calculating coefficients 

1. Light- cone gauge and non-renormalization theorem 

The calculations are most conveniently done in the light-cone gauge. This technology 
can be traced back to the pioneering work of 't Hooft ||, and has been well studied in the 
literature. In this formalism the energy-momentum vector is described by 

V ± = ±={E±p), (25) 
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2 

m . 



so that the mass-shell condition becomes p 2 = 2p + p_ 

Let us write down the Lagrangian of the model in the light-cone formalism: 



2i<9_ 



+ 9 j (E xlt a x)j (E x{t a x^j (26) 



In this formalism two-component quark g« fields are expressed via the one-component 
fermionic fields Xi, 




(27) 



(in the basis where 75 = 7°7 1 is diagonal). With the gauge fixed by A_ = 0, the A A 
component is expressed in terms of the quark fields. 
The weak interaction (|l0|) takes the form: 



rV 



G 



t f) a m i m Q 

— 



1 


t 






~~ r K~Xq 




jt XQ _ 



(28) 



A remarkable simplification occurs due to carrying light-like momentum q 2 = 
2g + g_ = 0. We can satisfy this condition by choosing the "spatial" component of the 
momentum q_ = 0, i.e., <9_0 = 0. (This means that the quantum is a left-mover.) Thus, 
on the <p "mass shell" the second term (containing <9_0) in Eq. (|28|) vanishes and the weak 
Qq<f) coupling takes a simple form, 



'"'weak 



G 

72n 



xIxq d+<p 



(29) 



Here we come to a very important point. In the 't Hooft model non-renormalization 
theorem for flavor non-diagonal currents at g_ = arises: all corrections to the weak 
vertex at this point vanish. In the light-cone formalism the kinematical point g_ = looks 
analogous to the zero recoil point in four- dimensional heavy quark theory. The analogy is a 
superficial one, however 0. In 1+3 QCD flavor non-diagonal currents at the zero recoil are 
renormalized by radiative and power corrections. 

To prove the theorem let us consider the gluon corrections to the weak coupling (p9"D- 
To order g 2 the relevant graphs are depicted in Fig. 0. 



In particular, for heavy-to-heavy transitions at q- = there is no dominance of the ground state 
production in the 't Hooft model. 
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FIG. 4. Radiative corrections to the weak coupling Qq<p. Dashed lines denote the gluons. 



The Feynman rules can be read off from Eqs. 
vertex is 



(25). In particular, the weak 



V = -v 



G 

/2n 



Q+ ■ 



(30) 



The graph f|a gives rise to following expression: 



AV = V ■ i 



Air 



d 2 k 



1 



2( Pq +k). 



+ k). 



(31) 



2(PQ+fc)- 



The integration over k + can easily be performed by the residue method. It is clear then 
that the integration over k + produces a nonzero result only in the case of opposite signs of 
(p q + and (pq + k)_ (the poles should be on the different sides of the integration path). 
For g_ = (pQ—p q )_ = these signs are certainly the same, (p q + k)_ = (pq + k)_, and there 
is no correction to the vertex. 

To finish up with the g 2 correction to the weak coupling we need to add graphs b and c 
of Fig. [| containing the self-energies £<g and S g of Q and q quarks. 



£q(p+,p-) 



Air J k 2 _ 



1 



p + + k 



(32) 



+ 2(p_+fc_) 



The integration contains a single pole only as a function of k + . Unlike the previously 
considered vertex correction, though, the integration over k + gives a nonzero result because 
the integral over the large semicircle in the complex plane of k + does not vanish. The 
integration over fc_ requires an infrared regularization. Following 't Hooft || we define 
the integration in Eq. (|32|) by putting a symmetric ultraviolet cutoff K for the k + integration, 
and a symmetric infrared cutoff A for the fc_ integration, 



\k\ > X . 



(33) 



4 OPE ensures that the dependence on the infrared regularization disappears in the width as long 
as it is the same at all stages of the calculation. 



13 



Then at \p + \ -C K the result for E is 



/3 2 /? 2 



2p_ 2A 



0(b-|-A). (34) 



The independence of S on p + means that no Z factor appears. The first term corresponds 
to a shift in the quark masses, 

m 2 Q ^m 2 Q -f3 2 , m 2 g ^m 2 q -f5 2 . (35) 

The second term produces a (noncovariant) shift in the reference point for the light-cone 
energy on mass shell, 

P + = ^— + 2A' (P - >A) - (36) 

This shift produces no effect on the widths. One-loop radiative corrections thus do not affect 
Qqcp transitions besides the mass shift given by Eq. (p5[). 

Moreover, it stays true for higher loops as well within the 't Hooft model. For in the limit 
iV c — > oo there are no fermion loop insertions into the gluon propagators. Then the higher 
loop corrections to the vertex, as well as to the self-energy, vanish in the way discussed above 
since the integration over k + yields zero. 

Notice that the non-renormalization theorem we derive within the 't Hooft model is a 
stronger statement than the one about zero recoil in four- dimensional QCD where radiative 
and power corrections break the non-renormalization of flavor non-diagonal currents. 



2. The leading coefficient cqq 

Now it is simple to account for higher orders in the coefficient cqq of the leading oper- 
ator QQ. To zeroth order in g 2 this coefficient was determined in Sec. ||, see Fig. |] and 
Eq. (Pp. As just discussed higher loop corrections merely shift the quark masses, Eq. (fto|, 
and therefore we get the coefficient cqq in all orders, 

G 2 171% — m 2 



Combining this result with Eq. (|20|) and with the suppression of the four-fermion operators, 



see Sec. |l 1 B 3| , we conclude that 



(i) there is no I/tuq corrections in the total width, much in the same way as in actual QCD 

[0,0; 

(ii) corrections l/nig, 1/jtoq and I/tuq to the total width are associated exclusively with 

the operator QQ. 
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In Sec. [Ty| we will prove a stronger statement: irrespective of the explicit form of these 
corrections, the hadronic saturation yields exactly the same result for the total width as the 
contribution of cqqQQ in OPE. 

The expression (|37| ) for cqq refers to a low normalization point, (3 <C // <C tuq. In order 
to calculate the matrix element of QQ over Hq we will need to express QQ in terms of 
Xq- It goes without saying that the operator QQ must be taken at the same normalization 
point. Then, the resulting series in j3 2 /trig cancels against a similar expansion coming from 
the operator QQ; there is no (3 dependence in the product cqqQQ. This can be seen by 
rewriting QQ in terms of the unrenormalized one-component field Xq an d mass uiq: 

QQ = Xq^Xq ■ (38) 

In evolving down to /i, higher orders lead to the substitution iiiq — > Jmq — (3 2 in this 
relation as well. With the quark mass substitution being the only effect of the radiative 
corrections we have 

Q2 j 

2 Im cqqQQ = —(m Q - w 2 ) XqJq~XQ ■ (39) 

The statement that the product cqqQQ is renormalization group invariant is trivial, of 
course. A nontrivial part of the result is encoded in Eq. (|39|) , which is valid to all orders in g. 
One could obtain this result by doing calculations at /i ^> rriQ when the mass of the Q quark 
coincides with its "bare" value rag, at fi = tuq, or at ji m<g, when a non-logarithmic 
evolution of the QQ operator and its coefficient functions must be taken into account, the 
outcome is the same, see Eq. (|39]). To make contact with the 't Hooft equation (i.e. to 
calculate QQ in terms of the 't Hooft wave function defined for bare quantities) we will need 
Eq. ( p9|) at the ultraviolet cutoff. Note that it can be conveniently rewritten as 

2lmcQ Q QQ = T QX 1 Q^XQ- (40) 



In Sec. [IV B| we will find the matrix element of XQ( m Q I '^-)Xq an d show that the corre- 



sponding expression for the total width coincides with the one obtained through the hadronic 
saturation. 



3. Four-fermion and multi-fermion operators 

As discussed above for the current-current weak interactions four-fermion operators do 
not arise to zeroth order in the strong coupling. Diagrams generating four-fermion operators 
in the g 2 order are shown in Fig. [|. 



15 



Q 



Q 



Q 




Q 



g 1 i 



a b 

FIG. 5. Four-fermion operators in g 2 order 

The light-cone gauge turns out to be again a convenient tool to use. Let us start for 
illustration with the simple one-loop diagram of Fig. || Feynman rules in the light-cone 
gauge, see Sec. |111 B 1| , lead to the following expression 



Imi J d 2 k5(2k + k. 



-k). 



(41) 



2(pq-&)- 



Here 5 (2fc + A;_) appears from the cut of the propagator, the cutting of the q quark propa- 
gator done by taking of imaginary part. The integration over A;_ is immediate due to k- = 
root of delta function (the other root, k + = 0, gives the same contribution, so just a factor 
2) 



/ 



Imi dk+ k-i 



1 



(PQ - k). 



71 



2(P0)- 



(Pq)+ ~ 



m„ 



2(P«)- 



(42) 



It means that the term 



const 



id L. 



2 1 



(43) 



appears in the transition operator. In the zeroth order in the 't Hooft coupling the equation 
of motion for the \Q ne ld is 



d+XQ = m 2 Q /(2id-)xQ ■ 



(44) 



In the rest frame zcL — » p_ = mg/i/2, and we reproduce Eq. ([TT|). 

Let us now apply the same technique to the loop part in Fig. [5|& in the limit of vanishing 
gluon momentum. Integrating over k_ we get 



Imi d/c_i_ k-i 



1 



7T . 



k)+- 



2(PQ)- 



It produces the term 



f i = const Xq^+Xq 



(45) 



(46) 



in the transition operator, note that the overall factor is the same as in Eq. fl43|). Summing 
up T and T% results in the substitution of id + by iD + = id + + A + . In this order the equation 
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of motion ( f44|) is also should be modified by the same substitution. Net result is that no 
change in OPE coefficients is produced by the loop in in Fig. [5j^. 

Note that it is true not only for vanishing momentum of gluon field but for any soft field 
as well, in other words, terms with derivatives of do not appear in the transition operator 
in one-loop order. Note also that if extra gluons are emitted from the loop we get zero for 
the diagram. Indeed, it increases the power in the integrand of Eq. fl45|) and the integral 
over the large semicirle in the complex plane of k + vanishes. It is the reason why diagrams 
of the type of Fig. |5]& give no rise to multi-fermion operators. 

To finish up with multi-fermion operators we need to account for diagrams of the type 
given in Fig. |5|a. The Feynman expression for this diagram after integrating over /c_ is 



Imi dk+ k 4 



1 



7T 



(Pq - k)+ - 



2(pq)- 



[(PQ ~ Pq) 



[(PQ-Pq)- 



(47) 



It is simple to check then that the diagram |5|a cancel out against similar diagrams where 
the gluon exchange is between q and Q quarks. In case of extra gluon insertions in one loop 
(relevant for six-fermion and higher dimension operators) we get a vanishing result right 
away. 

Thus, we proved that four-fermion and multi-fermion operators do not arise at the level 
of one loop. They show up at level of the second loop, see the diagram [| for four-fermion 
operator. The dimensional counting reveals then that four-fermion operators give l/m^ 
correction to the total width. 











g' 1 


'g 


q. : i 


: qi 


FIG. 6. Four-fermion operators in g 4 order 



C. OPE representation for inclusive width 

Putting everything together we get the OPE representation for the inclusive width: 

( H q\Xqj^Xq\Hq) 



G 2 mi 



Air 



rn 



Q 



(Hq\QQ\H q ) 



2M, 



2M 



O 



Hq 



We have thus obtained a very simple result: 

• The partonic expression Tq represents the asymptotic term for itlq — > oo. 



m 



Q. 



(48) 
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There is no I/tuq contribution in OPE as long as the weak interactions are of the 
V x V type. 

Through order l/rrin only a single operator contributes, QQ. 

The leading correction ~ 0(l/mq) enters through the expectation value 
(l/M HQ )(H Q \xUm Q /id-)xQ\H Q ). 



The second part of Eq. (IS) has been written in terms of the light-cone operators to provide 



a way of rewriting the matrix element in terms of the 't Hooft wave function of the hadron 
Hq. The OPE result for the inclusive width can be recast in terms of the sum over exclusive 
hadronic channels. This will be proven next. 



IV. MATCH BETWEEN OPE-BASED EXPRESSIONS AND HADRONIC 

SATURATION 



A. Exclusive widths via the 't Hooft wave function 

With 1+1 QCD describing manifestly confining dynamics, its spectrum consists of 
mesonic quark-antiquark bound states. In N c — > oo limit these mesons are stable in re- 
gard to strong decays. The masses and the light-cone wave functions (p(x) (with xe [0, 1] 
having a meaning of a portion of momentum carried by the quark) of these mesons can be 
determined as eigenf unctions and eigenvalues of the 't Hooft equation. In particular, the 
initial state Hq = [Qq sp \ is the ground state in the sector with the the heavy quark Q and 
the spectator antiquark q sp . Its wave function ip Hg satisfies the following equation: 



M 2 Hq <Ph q (.x) 



m, 



P 2 



x 



+ 



sp 



P 2 



1 — X 



Jo (y - xy 



(49) 



where m sp denotes the mass of the spectator antiquark and the integral is understood in 
the principal value prescription. The solutions to the equation are singular at x = and 
x = 1 where their behavior is given by x 70 and (1 — x 
the following conditions: 



_JHo_ 

tan7T7o 



m 



Q 



P 2 



2 



tan7T7i 



71 , respectively, with 7 01 defined by 

P 2 



m 



sp 



P 2 



(50) 



The masses M n and wave functions (p n of final mesons h r 
't Hooft equation with uiq substituted by m q : 



M 2 M%) 



P 2 



+ 



m 



sp 



P 2 



x 1 — X 

The functions <p n form a complete basis, i.e., 



[<l<lsp\n are defined by the same 

<Pn(y) 



<p n (x) - P 2 dy 

Jo {y — x) z 



^2<p n (x)(p n (y) = 5(x - y) 



(51) 



(52) 
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Let us find out now how an exclusive width T n of Hq — > /i n decay is expressed via wave 
functions <Ph q , fn of initial and final mesons, 



2M Hq {Ml Q -MD 



G 

(Ml xjxg 09+ 1 #q ) 



M£ Q - M n 2 
4tt M# 



(K\x\xq\h q ) 



(53) 



where the factor l/(Mfj Q — M 2 ) is the Lorentz invariant phase space (LIPS) of two-particle 
final state and for the matrix element we have used Eq. (|29|) in the kinematics where the 
momentum is 



q- = 0, 



1 



V2M, 



H, 



M 2 Hq - Ml 



(54) 



It is simple then to write the matrix element of xIxq m terms of the 't Hooft wave functions, 



G 2 Ml Q - Ml 
4tt M Ho 



dxip n (x)(p HQ (x] 



(55) 



B. Sum rules 



Using the completeness condition (J52|) we can derive sum rules for these partial widths 
by weighing them with powers of Mfj — Ml. The first one is 



AtcM 



E 



L, 



E 



G 2 foM 2 HQ -Ml n=0 



&X(p n (x)(p HQ {x) 



(56) 



Note that the sum runs over all states h n including those unaccessible in the real decays of 
Hq, i.e. with masses M n > Mh q - These transitions are still measurable by the process of 
inelastic lepton scattering off Hq meson. This sum rule is an analog of first Bjorken sum 



rule and was discussed in (18| . 

To get next sum rules let us multiply Eq. ( f49"D by <p n (x) and Eq. (|5lD by <p HQ (x), respec- 
tively. After integrating over x and subtracting we find: 



(M 2 Hq - Ml) / dx ip n {x)(f HQ {x) = {m 2 Q - m 2 ) / —<p n (x)ip HQ (x) . 

«/ J X 



(57) 



Two more sum rules then arise: 



n=0 



G 2 



\m Q - m qJ 



2 , / 1 dx o 



X 



(58) 



4vrM 



He 



G 2 



oo 

E 

n=0 



r„(M 2 



Hq 



Mi 



m 



Q 



(59) 
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The second and third sum rules differ from the first one in two aspects: they depend on 
the quark masses mq and m q explicitly and the integral over the wave function is not fixed 
by a normalization condition; it can, however be calculated in the 't Hooft model. One 
should note that while the integrand (p 2 H (x)/x 2 ~ x _/3 ^ m «, is singular at x = 0, it is still 
integrable, since [3 2 /mq <C 1. Note also that expanding the second sum rule ( |58D in 1/mq 



to the linear order we reproduce the corresponding sum rule of Ref. [IH|. 

The sum rules above provide us with detailed information on the saturation of the sums 
over the final hadronic state. The quantity 

Wn = G 2 Ml -Ml (60) 

can be interpreted as a normalized probability of producing the state n. Indeed J2 w n — 1 
according to the first sum rule (J56D - Then the sum rule (|58D implies 

(M| Q -M^) = K-^)/ o 1 ^^ Q (x). (61) 

If both masses m sp and m q are smaller or of order of (3 we conclude that 

Ml) = ^-l^ m 2 Q + 0{(3 2 ) ~ (3m Q . (62) 

Here we have anticipated the result for (l/x) from Eq. ([71]) in the next subsection, in 
conjunction with Eq. ([7(]). The reason for (M 2 ) ~ (3mq is clear on the physical grounds: 
in the partonic approximation the final state is formed by the quark q with the momentum 
mq/2 and by the spectator antiquark q s p with the momentum of order f3. 

The sum rule (^) (after subtracting the square of the second sum rule (]55|)) determines 
the dispersion: 

(Mt) - (Ml)' = (1) - (i) 2 ~ {Ml)' . (63) 

What about higher moments? It is not difficult to see that the next one, Y^Lo ^n{Ml — 
M 2 ) 2 , is a divergent sum because l p 2 H n {x)/x' i would no longer be integrable. It defines the 
asymptotics of T n at large n, 

r.ocjL. (64) 

n 

Let us recall that M 2 = TT 2 (3 2 n for high excitations, n 3> 1. We will see in the next 
subsection that the asymptotics (]64) matches the contribution of the four-fermion operators 
in OPE. 



C. Matching 



Armed with exact sum rules (pq) , flog), (p9|) we are well prepared to verify a perfect 
match between the OPE result and the expression obtained by summing over hadronic final 
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states. From the second sum rule (pg) we have for the total semileptonic width Th q the 
following relation 

rB Q = ¥' m \7 mi f E r «- (65) 

4 ^ Hq Jo 1 Mn ^ HQ 

The second term is actually positive (r„ defined by Eq. d53|) is negative at M n > M Hq ). 
Using Eq. (|64D its size can be estimated 

- r « K — T ( 66 ) 

M n >M„ Q m Q 

Thus we have derived from the 't Hooft equation the following result for the inclusive 
width 



^ [1 -V 2 Ho {x) + 0\^-\ . (07) 




This expression coincides with the OPE result of Eq. fl4q ) as seen by rewriting the matrix 
element in Eq. ([48]) in terms of the ground state wave function </?# Q (x), 



( H q\Xq^Xq\Hq) m Q f 1 dx 2 

■¥>h ( x ) ( 68 ) 



2M Hq M Hq Jo x ru i 

where, besides normalization factors, we have used also the substitution id- — > xMh q / V%- 
This completes the proof of the perfect matching between OPE and the hadronic sat- 
uration through the order 1/rriQ. Let us stress that the matrix element fl68D given by the 
integral over the ground state wave function is implicitly uiq dependent, its leading term 
is 1 followed by 1/rriQ and higher terms (see the discussion in Sec. |III A| ). In the 't Hooft 
model one can, of course, evaluate the matrix element explicitly although it is not relevant 
to our main objective - probing the quark-hadron duality. 

The absence of 1/ tuq corrections was demonstrated in Sec. [IH A| by operator methods. 



Let us show now that the same statement can be derived from the 't Hooft equation as well. 



To this end we use the approach of Ref. Jl8| to the heavy quark limit, generalizing it to 



include l/mn corrections. Instead of x, the appropriate variable for the large mass limit is 

t=(l-x)m Q , ip(x) = y/mQ(t)(t) (69) 

An expansion of the 't Hooft equation in I/trq (after a substitution of the new variable) 
and the virial theorem lead to the following relation for the meson mass, 

^ = 1 + 2 «_4 +0 (^), (70) 



we 



where averaging is over the Hq meson wave function 4>(t). As compared with Ref. [K 
have added a (3 2 /rriq term. Its origin is simple: it accounts for the renormalization of the 
heavy quark mass. Note that (t) ~ (3 provided m sp < (3. 
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We also need a similar expansion for the integral entering Eq. (|67j), 



'^ W = 1+ W + W +0 (M. (71) 

o x o m Q m 2 V m o/ 



Substituting Eqs. ([ft], [TT|) in Eq. ($7$) we have, 

l + i(/3 2 -(t 2 ) + (t) 2 ) -COI^-I . (72) 



- i 




This expansion should be compared with the operator representation (PU[). The I/tuq ex- 
pansion of cqqQQ produces the same 1/m.g term; the part oc /3 2 comes from the expansion 
of Cq Q , see Eq. fl37[). 

Note that without the /3 2 term the correction to 1 is negative, i.e. 



2m 



K 7^ • (73) 



One more comment about 1/rng terms. In the OPE approach they are due to the 
four-fermion operators generated by the graph in Fig. ^. Although the corresponding OPE 
coefficients are not calculated, the consideration above shows that the contribution of the 
four-fermion operators is dual to the sum of T n with M n > M^ Q for final state mesons, i.e. 
channels kinematically inaccessible in the decay. 



V. VIOLATIONS OF DUALITY 

A. Global and local duality 

Having established a perfect match between the OPE prediction for the total width and 
the result of the saturation by exclusive decay modes, through 0(1/i71q), we must now turn 
to the issue of where the OPE-based prediction is supposed to fail. The failure usually 
goes under the name of "duality violations", a topic under intense scrutiny in the current 
literature. The definition of what duality violation is varies from publication to publication. 
Quite often, the researchers in the field stick to a vague notion of deviations "of certain rates 
for processes involving hadrons from the underlying partonic rates" . This is, for instance, the 
convention of Ref. |10| where duality is understood as the coincidence with the parton-model 



prediction. If so, any nonperturbative contribution to the given rate would be interpreted 
as a "duality violation" , which does not make much sense to us. 

We must precisely define what is meant by duality and its violations. Assume that a 
certain process is amenable to calculations within OPE. This means that an appropriate Eu- 
clidean quantity can be chosen, and the OPE series can be constructed. This series presents 
the quantity of interest as an expansion in an inverse large parameter, e.g. 1/Q 2 or 1/E. The 
very same quantity can be expressed as a dispersion integral over the imaginary part defined 
in Minkowski space. In e + e~ annihilation the imaginary part coincides with R(e + e~), in the 
transition amplitudes for heavy flavors the imaginary part reduces to semileptonic spectral 
densities, etc. 
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In order to treat the nonleptonic decays in the same vein one can introduce a spurion 
in the weak vertex, carrying a momentum q. In other words, let us substitute the weak 
lagrangian by 



where S{x) is a spurion field. Now let us consider the forward amplitude A(s) of the process 



as a function of s = (p + q) 2 . This variable s plays the same role as s = q 2 in e + e 
annihilation. The total nonleptonic width is given by Im^4 at s — MjL. We are free 
to consider Im^4(s) in the complex s plane where it has two cuts: at s > and at s < 
—2Mjj Q + 2q 2 (the second cut is due to u channel). Choosing a reference point so far away 
from both cuts but closer to the first one we can express A(so) as a dispersion integral over 
the discontinuity across the cuts. On the other hand at the very same point so one can apply 
the operator product expansion for calculating A(sq) in terms of matrix elements of local 
operators, (HqIO^Hq) . This gives sum rules which allow us to determine Im^4(s) at large 
s, in particular, at s — Mjj Q . This is fully analogous to what one does in e + e~ annihilation 
for R(s). In both cases smoothness is assumed (of course, in in e + e~ annihilation all positive 
values of s are accessible and one can check this assumption while in the case of nonleptonic 
decays Mh q is fixed). Needless to say that the total semileptonic widths can be treated 
along the same lines, the only difference is the presence of leptons in the intermediate states. 

In the semileptonic decays OPE allows to predict, additionally, various distributions in 
the lepton momenta. This is, probably, the reason why it is usually claimed that the status 
of duality is more solid in the semileptonic decays. To show that it is not the case let us 
consider the semileptonic decays with the light quark in the final state. The OPE-based 
predictions for the spectral distributions are valid almost everywhere; they fail only in the 
end-point domain Jl5|| . For this reason the total semileptonic widths cannot be obtained 
by integrating over the spectrum if we want a prediction which includes the linear in I/tuq 
corrections, it is why the argumentation in Ref. |H| was not sufficient to prove the absence 
of 1/rriQ corrections. Nevertheless, the statement of absence of such corrections in the total 
semileptonic widths can certainly be justified by the procedure described above for the 
nonleptonic widths. Thus, we see that the theoretical status of all these processes - e + e~ 
annihilation, semileptonic and nonleptonic decays of heavy flavors - is basically the same. 

By performing an appropriate expansion of the dispersion integral we obtain sum rules 
relating certain moments of the imaginary part of transition amplitude to matrix elements of 
consecutive terms in the OPE series constructed in the Euclidean domain. The predictions 
obtained in this way will be referred to as global duality. Taken at their face value, they are 
exact, to the extent we can calculate the coefficient functions and matrix elements of the 
operators involved in OPE. No additional assumptions are made. The predictions obtained 
in this way are consequences of fundamental QCD. Therefore, it does not make any sense 
to speak about violations of the global duality. One can only speak of the precision of 
calculation of the coefficient functions and determination of the matrix elements. 

Unfortunately the term global duality is often used in a loose and ambiguous sense. It 
is applied indiscriminantly to integrals over the spectral densities with the weight functions 




(74) 



S(q) + H Q (p) — > light hadrons — ► S(q) + H Q (p) 



(75) 
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chosen ad hoc. Our definition is narrower: it refers only to those specific integrals which 
emerge from the dispersion representation. 

The notion of local duality on the other hand requires further assumptions. Assume that 
we want to predict imaginary parts (spectral densities) point by point, at large energies (or 
q 2 ). If one assumes that the spectral densities at the given energy are smooth, then from 
the moment integrals we can certainly predict the densities themselves. This amounts to an 
analytic continuation of the OPE series (truncated in a certain way), term by term from the 
Euclidean to Minkowski domain, with the subsequent calculation of the imaginary parts of 
each individual term in the series. The prediction obtained in this way is evidently a smooth 
function of the parameters. We then compare this prediction with the quantity measured in 
terms of hadronic contributions. The difference between the OPE-based smooth result and 
the experimental hadronic measurement is referred to as the duality violation meaning the 
violation of local duality. 

By its nature the OPE results are series in powers of A QCD /_E and do not account for 
terms like exp[— (_E/Aq CD ) fc ] (in the Euclidean domain). Although such terms are due to 
large distances, a signal of their appearance could show up in the short distance OPE series 
in the form of a factorial divergence of the series in higher dimensions. The situation is 
reminiscent of that in the perturbative expansion. The divergent a s series (e.g. due to 
infrared renormalons) give rise to terms exp(— C/a s ) although such terms can appear even 
in the absence of renormalons (for instance, as the quark condensate). 

In other words, the OPE construction accounts properly for short distance singulari- 
ties while the exponential terms are due to large distances being nonsingular at short dis- 
tances. Thus, the duality violation is something we do not see in the (truncated) OPE series. 
The duality-violating terms are exponential in the Euclidean domain and oscillating (like 
sin[(_E/AQ CD ) fc ]) in the Minkowski domain. 

From this standpoint there is no distinction between, say, the total e + e~ annihilation 
cross section or the semileptonic rates of heavy flavors, on the one hand, and the nonleptonic 
rates of heavy flavors, on the other. Sometimes it is claimed that the former processes are 
"pure" while the latter are "impure"; it is even asserted that "duality follows from OPE 
in the first case while it has no theoretical justification in the second case". We assert 
that a duality violating exponential/oscillating component, associated with the neglected 
tails of OPE, is present in all processes, and the only physically meaningful question is its 
magnitude, as a function of large parameters (e.g. a momentum transfer or toq) and specific 
details of the process under consideration. 

B. Oscillating terms in 't Hooft model 

The appearance of duality violations in the form of oscillating terms is evident in the 
't Hooft model where the spectral density is formed by zero-width discrete states. Indeed, 
each time a new decay channel opens dY n Q j drriQ experiences a jump (Th q is continuous) , so 
that immediately above threshold dT h q / drriQ is larger than the smooth OPE curve, in the 
middle between two successive thresholds it crosses the smooth prediction, and immediately 
below the next threshold dT HQ /dmQ is lower than the OPE-based expectation. 
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The amplitude of oscillations can be estimated as follows. Let us present the total width 
T Hq as 



n=0 M n >M Hn 



(76) 



Widths T n are exclusive widths of two-body decays, Hq — » <fi + h n , where h n is the n-th 
excited qq sp state with the mass M n . For M n > M Hq widths Y n are not, of course, physical 
ones for Hq decays but they are well defined. We have used this presentation in Sec. [IV (J| 
where it was shown that the first term can be related to the wave function of Hq state on 
one hand and to the matrix element of QQ operator on the other. It is clear that the first 
term in Eq. (|76j) is a smooth function of vhq and contains no nonanalytic terms we are going 
after, they are in the second one. 

Thus we need to to know T„ for M n in the vicinity of Mh q . For M n 3> Mu Q we found in 
Sec. |IV B| that T n oc M~ 6 , see Eq. fl64|). To extrapolate to the vicinity of M Hq we account 
for the threshold factor, T n oc (Mfj Q — M%)/M^. The coefficient in this dependence can be 
fixed by duality of the second term in Eq. (|76|) to four-fermion operators 0^ q , discussed in 
Sec. WQ , 



M n >M Hn 



[H Q \O iq \H Q ) 



2M, 



He 



(77) 



We did not calculate coefficients but we found that up to numerical factor 

G 2 (3 A 
2n rrin 



(78) 



Assuming m sp , m q ^ (3 we estimate (HqIO^Hq) /2Mh q ~ (3. Then our model for r„ in 
the range M n > Mh q is 



n M 2 H -Ml 

y M 8 



(79) 



The sum in in Eq. ( |77D was approximated by the integral using M 2 = ir 2 f3 2 n. 

Now we can evaluate the sum (|77|) with the better accuracy accounting for small nonan- 
alytical terms. The result is 



- ^ r n = Ar OPE + Ar° 

M n >M Hn 



G 2 (3 
2vr 



P 



7T 4 / (3 



M 



Ar° sc = ^-tt 3 G 2 (3 



HqJ 



2 M. 



M. 



Hq, 



x(l - x) - - 



(80) 



where 
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x = Fractional Part of —£§r , x G [0, 1) . (81) 

y 7i 2 p 2 j 

The smooth part Ar OPE in Eq. (]80f) is given by same OPE term of Eq. ( jfTj) we discussed 
above plus higher in the power of 1/ttlq corrections . The nonanalytic in Mfj Q part Ar° sc 
oscillates with period vr 2 /3 2 , see its plot in Fig. [7]. 
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FIG. 7. Oscillations in T^ Q . The ratio AF° SC /G' Z P is presented as a function of Mj^/vr 2 /? 2 



The amplitude of oscillation is 



Ar° 



3tt 4 







M. 



Hq, 



(82) 



Note that the derivative d(AT° sc ) / drriQ contains discontinuities at thresholds, the amplitude 
of oscillations is larger for the derivative, 



d(Ary sc /dm Q 



dT Q /dm Q 



12n 2 







M. 



(83) 



The oscillations under discussion cannot be produced by any truncated OPE series, they 
are not seen in the OPE. Thus the estimate fl8"2|) gives the actual scale of the expected duality 
violations in the problem at hand. Of course, this estimate is obtained within a specific model 
for r n . The gross features of the result are independent of the model, however. They are 
determined only by the fact that the resonances have zero widths. 

It should be noted that the limit iV c — > oo presents a scenario which maximizes duality 
violations. In this limit the thresholds open "abruptly", right at the position of the reso- 
nances, since the resonance widths vanish. In the real world of finite N c the highly excited 
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states have finite widths, and this effect, on its own, smears the hadron-saturated cross 
sections dynamically. If in the zero width approximation the oscillating duality violating 
component is suppressed only by powers of a large parameter (1/rriQ in the case at hand), 
switching on finite widths will further suppress the oscillating component exponentially, see 
Sec. IV D| and Ref. 



12 for further details. 



C. Lessons 

From the considerations above we conclude that: 

(i) Oscillating terms violating local duality are definitely present in the total decay rate 
(considered as a function of tuq); 

(ii) Amplitude of oscillations is ~ 0(l/mV), i.e. strongly suppressed; 

(iii) If we could average over rug in a sufficiently large interval the power suppression of the 
oscillations would turn into exponential suppression (in actual QCD, with N c = 3, the 
finite resonance widths do a similar job). Then it is perfectly legitimate to consider 
the OPE-based predictions beyond 1/mg. 

With this understanding in mind we now turn to the discussion of duality violations in 
actual four-dimensional QCD. 



D. r decays in 1+3 dimensions 



Let us discuss a quantity of practical interest in 1+3 dimensions along similar lines, 
namely the normalized hadronic r width R T : 

R T = rCr-^ + hMirons) 
r(r -> v T e-v e ) 

It can be expressed in terms of spectral densities py and pa in the vector and axial-vector 
channels, respectively, 

J 2 (M 2 ) , J 3 (M 2 ) 



P [ M ? ds ( sYf s\ I (M 2 ) I 2 (M 2 ) 



(85) 



where the moments I n are defined as 



r M 2 

I n (M) = dss n [ Pv (s) + p A (s)] . (86) 

Jo 

While t decays represent a simpler dynamical problem than the weak decays of heavy flavor 
hadrons we have to simplify it further still before we can arrive at some definite conclusions. 
To estimate the oscillating contribution to R T which constitutes duality violation that cannot 
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be seen in a truncated OPE we consider the limiting cases of M T and N c large. We will 
show that, adopting the resonance model motivated by two-dimensional QCD, for N c — > oo 
and M T large, yet finite, the duality violation in R T scales as 1/M®; for M T — > oo with N c 
large, though finite, the oscillating term is suppressed exponentially. 

Our consideration will be admittedly illustrative. One should not take literally the 
numbers we will obtain for many reasons: first of all the r mass is not much larger than the 
spacing between the resonances, second, N c is not large enough to warrant the zero width 
approximation. Still we believe that the consideration is instructive in a qualitative aspect. 

For large N c the spectrum of 1+3 QCD is expected to consist of an infinite comb of 



narrow resonances - in complete analogy to the 't Hooft model [2U|. To keep the closest 
parallel to it we further assume that the high excitations in a given channel (e.g. the vector 
channel) are equally spaced in m 2 . This agrees with the general expectation of a string-like 
realization of confinement leading to asymptotically linear Regge trajectories. The masses 
of the excited states in, say, the p channel are then givenf] by m 2 = m 2 + 2n/a' PTf , 
with a' being the slope of the Regge trajectory (for a review see p2|). Experimentally one 
finds 2/o/ ~ 2 GeV 2 . For large values of s the spectral densities for both the vector and 
axial-vector channels will approach the form: 

Pv (s)=p A (s)=N c .j:s(--n); a 2 = -, (87) 



n=l ^ / « 



where a special notation a 2 is introduced for 2/ a!. Equation (K7|) is clearly not expected 
to hold at moderate and small values of s where the vector and axial-vector channels are 
drastically different and the resonances are not equidistant. However, details of the spectral 
densities at small s play no role in duality violation. 

The contribution of any particular resonance of mass M/, to R T , according to Eq. (jS5|), 
is given by a simple polynomial in 1/M 2 (times the step function Q(M 2 — M|) ). Therefore, 
variations of properties of resonances below a certain fixed mass p change only the regular 
terms of the 1/M 2 expansion, but have no impact on the oscillatory component. From 
Eq. fl8~5|) it is clear that such variations change only coefficients of the 1/M 2 , l/Mj 3 and 
1/M^ terms. It will be clear in what follows that the formal OPE for R T exactly reproduces 
these three expansion coefficients as well. 

The spectral density in Eq. ( |S7|) is dual to the parton model result; i.e., it coincides with 
it after averaging over energy, 

(p VjA (s)) = N c . (88) 

Thus, the asymptotic prediction for R T at M 2 — > 00 is 

R° T = N c . (89) 

The sum over resonances in R r is easily calculated analytically: for the spectral density 
of Eq. flgTj) it is 



5 In other, less QCD-friendly scenarios, one obtains instead m 2 = m 2 + n/a' . The distinctions 
between these two scenarios are irrelevant for our discussion. 



28 



where 



R T 
R OPE 

IvT 

5R osc 



R°™ + SR osc 
- 2 1 



1 - 



a 



a- 



Ml 30 \ Ml 



-x(l - x)(l - 2x) 



a 



2 \ 3 



Ml 



+ 



x 2 (l-x) 2 - 



30 



Ml 



x = Fractional Part of — j- , x E [0, 1 



Ml 



(90) 



We presented the result as a sum of two functions of M 2 , the first one, R° PE , is a smooth 
function expandable in 1/M 2 . The second one, 5R osc , oscillates with the period a 2 ; its 
average vanishes, see the plot of 5R osc /R® in Fig. |^. 
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FIG. 8. Oscillations in R T . The plot of 5R osc /R® is presented as a function of M 2 / a 



Let us show now that i?° PE coincides with the OPE prediction in the model. Power 
corrections can be presented in the following way: 



(91) 



where the "condensates" /„ are 



dss n Ipv(s) +p A (s) -2N c 



(92) 



This integral representations for the "condensates" I n follows from Eqs. (|S5D, (|86| ) if one 
assumes that the spectral densities approach their asymptotic limits faster than any power 
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of 1/s. In the model at hand, with the comb-like spectral density, the integral representa- 
tion (^) requires regularization. As a regularization one can introduce the weight factor 
exp(— es), taking the limit e — > at the end. With this regularization, i?° PE from Eq. fl9~0"|) 
is reproduced. 

To justify the procedure it is instructive to make one step back and follow more liter- 
ally the original OPE procedure. Namely, our primary object of interest is the Euclidean 
polarization operator 



n(Q 2 



7T 



ds 



pv{s)+ p A {s) 
s + Q 2 



(93) 



The original meaning of I n is the coefficients in the asymptotic expansion of II (Q 2 ) at 
Q 2 -> oo : 



7rn(Q 2 ) = - 
The comb of model of Eq. 



2N C In Q 2 + const + £ (-l) n 

n=0 



{Q 2 ) n+1 



(94) 



was considered in Ref. 



vrn(g 5 



1 



\ (Q 2 /cr 2 ) + n 



+ const 



-2N C 



Q 



+ const 



(95) 



where if) is Euler's if) function. The asymptotic expansion of H(Q 2 ) in the model takes the 
form 



vrn(Q 2 



-2N P . 



91 \± _ v ^ 

11 v 2 2Q 2 ^ 2n \Q 2 / 



+ const , 



(96) 



where B n stand for the Bernoulli numbers. Taking the coefficients of 1/Q 2 , 1/Q 6 and 1/Q 8 
terms from this equation we find the consistency with i?° PE of Eq. (0). (The term 1/Q 6 is 
absent and 1/Q 8 is defined by B A = —1/30.) 

Note that the spectral density ( jB7f) is not literally "QCD-compatible" : the corresponding 
n(<5 2 ) has the 1/Q 2 nonperturbative correction forbidden in QCD II. This can be easily 
cured by adding the resonance with n = with a half weight, which just would amount to 
adding N c a 2 /Q 2 to nIl(Q 2 ). Since this obviously does not change SR osc at any value of M T , 
this is inessential for us. 

Let us discuss now the duality- violating 5R osc (M 2 ) (see Fig. |^). Its dominant compo- 
nent scales as 1/Mf . It is intriguing to note that the very same scaling law was obtained in 



Ref. |23[ from totally different considerations invoking instantons. This oscillating compo- 
nent vanishes at M T corresponding to the new thresholds, and at one point in the middle 
between the successive resonances; only the second derivative of R T has a jump at the 
thresholds, 



1 



d 



x=0 



R T IdM 2 



R T 



(97) 



x=l 
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This is the consequence of the threshold factor (1 — s/M 2 ) 2 in Eq. (|85|). One power in it is 
just the two-body phase space factor \p\/M T . In 1 + 1 dimensions one would have instead 
l/|p|, the threshold singularity, and the duality- violating component would be enhanced 
correspondingly. 

The average of 5R osc vanishes while the amplitude of oscillations amounts to 



5R osc 



R T 



2 \ 3 



a" \ 



- 



3VT2 \M 2 ) 



(98) 



It is interesting that, in spite of the fact that 5R osc is given by a polynomial between the 
thresholds, the whole function is well approximated by 

1 / M 2 \ fa 2 \ 3 

2tt-^ t« • (99) 



3V12 

Note that the numerical coefficient in Eq. fl9lp is rather small, compare e.g. to Eq, 
This suppression is related to the fact that the characteristic scale is o" 2 /27r rather then a 2 . 

Taking our estimate of the oscillation amplitude at its face value and using the actual 
value of the r mass in Eq. (|98|) we find 5R osc /R T ~ 3%. It is clear that this is a very 
crude estimate given the fact that in the actual r decay a 2 /M 2 ~ 2/3 and we deal with one 
oscillation at most. 

In the real world with N c = 3 we expect a further suppression of deviations from duality 
due to the nonvanishing widths of the resonances naturally smearing out the amplitude of 
the oscillation. A rough estimate of this effect can be given in close analogy to Ref. Jl2| . 
Let us introduce a dimensionless constant B representing the width-to-mass ratio: 

^ = | (1 + 0(1/A y ); (10 o) 



i.e., B stays finite for large N c . One actually guestimates B ~ 0.5. Then we infer (see [12 
for details) 

AR T 1 fa 2 V ( 2nBM 2 \ , s 

The power-suppressed oscillations eventually turn into exponentially suppressed, although 
at a larger energy. 



VI. NONVANISHING AND COMMENTS ON THE LITERATURE 

The work [[K|] stimulated our interest in the 't Hooft model as a laboratory for exploring 
heavy quark expansions in inclusive decays, and the implementation of duality. The authors 
of Ref. ]10j compared the decay width of a heavy flavor hadron in the parton approximation 
with the result obtained by summing over all exclusive transition rates for Hq — > hjij. A 
systematic excess of the total width Y Hq over its parton value Yq was observed and was 
fitted to be 
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T Hq -T q 0.15 



(102) 



The authors interpreted this violation of duality. According to our understanding 

(see Sec. |V A|) it should rather be called breaking of OPE. 

Our analytical treatment does not support such a conclusion. In particular, we have 
proved the absence of linear in 1/mq corrections not only by the OPE method but by the 
direct analysis of the 't Hooft equation (see Sec. |TV| ). Although we tried to closely follow 
the analysis of Ref. | [10|] , still there is a difference of kinematical nature. We considered the 



ip fields (leptons or quarks) to be massless, while = m q ^ in Ref. ||10|| . The choice of 
m^p = allows us to limit ourselves to the point q 2 = where great simplifications occur, in 
particular, only the massless pion state is produced in the ipip channel. Analytic solution of 
the problem turns out to be possible. 

We have checked that the analytical expression for the triple overlap integral of Ref. [K^ 



reduces at q 2 = to a simple overlap integral ([55]) of two wave functions. The simple 
structure of (|55"D combined with the completeness of the wave functions is sufficient to prove 
a perfect match between OPE and hadronic saturation, at the level of l/trig. 



To make full contact between our results and those of Ref. [T(| we must now consider 



the impact of ^ 0. Needless to say that at ^ the parton result for Tq changes, 



T Q (m^ ^ 0) = T Q (m^ = 0) 



(103) 



This effect was certainly included in the analysis of Ref. | lQfl . It is also accounted for in 
OPE as a change of the coefficient of the operator QQ. 

The leading effect due to ^ is linear in m^, however. Thus, one can wonder whether 
it produces 1/rriQ corrections. In the next subsection we will show that the linear in 
corrections to the total width are suppressed as l/rrVq. Moreover, this is a leading effect 
which produces a distinction between semileptonic and nonleptonic total widths, all other 
effects which differentiate them are suppressed as I/itlq. It is clear then that the 1/rriQ 
violation of duality claimed in Ref. |l(| if it would be present in the nonleptonic width must 



have been present in the semileptonic width as well. 

Corrections I/tuq come also from quadratic in terms in OPE. As it was mentioned 
above these corrections do not differentiate semileptonic and nonleptonic widths. We esti- 
mate them in Sec. VI B . We also estimate in Sec. [VI Q effects of novanishing for violations 



of local duality. Overall, we conclude that the OPE approach shows that the nonvanishing 
results in 1/ m?q corrections which are numerically small and cannot explain the alleged 
deviation (11021). 



On the hadronic saturation side we have checked that the triple overlap integral of 
Ref. [plj is expandable in m^, and the leading correction is quadratic in m^ >q . Assuming 
that m^p < (3 and performing an analytic summation of the widths using the expressions of 
Ref. [[U| for the amplitudes in conjunction with the sum rules for the hadronic polarization 
tensor Imll(g 2 ) similar to those in Eqs. (|^,^,^), we obtained exactly the same m^/mg 



correction as in Eq. (|103|) . This holds for the inclusive width smeared over a small interval 
of masses as discussed in Sec. |V1 CJ Moreover, we checked the exact matching with OPE at 
the level 1/itiq. 
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A. Linear in mv, corrections 



All effects due to 7^ reside in 11^^, see Eq. (|7|). The linear in part in OPE for 
IV is 

An^( ff ) = ^ (o|^|o) (Ml _ ^ . ( i 04 ) 

It is easy to see that this effect is nothing but a shift of the pion (0 = [V'V'Jo) mass from 
zero, 

f4 = - W m f (0\W\0) . (105) 

Indeed, comparing Eq. @ and Eq. ( |104| ) we see that the latter is the first order term in /x^ 
expansion of the pion pole l/(q 2 - fify in = -(l/7r)(g M g v - q 2 g iilV )l (q 2 - fity. Thus, in 
order to take into account linear in effects in the transition operator all one needs to do 
is to calculate the decay Q — > <ft + q with the nonvanishing pion mass. The result is 

ATfi 2mnm a 8ir . . - , . , mnm a 

T Q ( m Q - m q) N c \ m Q ~ m q) 

At large tuq it falls off as I/tuq. An extra suppression m q /mQ is specific for two dimensions. 

As for for the numerical value of the correction ( |106|) we we will substitute {4>ip) = 
— jV c /?/\/l2 in the chiral limit [HJTQ]. For the value of = m q = 0.56/? adopted in 



Ref. |10j we get 

Ar (i) 

—^ = 2.271 — 1 . (107) 




B. Quadratic in r/u corrections 

The analysis of the previous subsection refers to the operator QQipif) in OPE for the 
total width. Another operator generated due to ^ is the four-fermion operator of the 
type QQqq. It appears from the graph of Fig. [5] with lepton lines substituted by ip lines. 
Calculation of this graph is a simple exercise. The result for the corresponding correction 
to the width can be presented in the following form 

Ar (2) G 2 ml {Hq I Qj^q q~if,Q \ Hq) 

Ar «« - N ^ Q mr Q — ■ (108) 

Up to the factor N c the same contribution appears in the semileptonic width. In difference 
with the operator QQipip the large N c limit does not allow us to factorize the matrix element 
of Eq. (|108|) . We will use the factorization to get an estimate, 
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[H Q \Qi 5q(n5 Q\H Q } ~ -JL(0|gg|0) (H Q \QQ\H Q ) . (109) 



As in the previous subsection we use the chiral limit value (qq) = — N c f3/yV2 of the 
TT|JT9|| and = m g = 0.56/3 for the numerical estimate, 



vcv 



0.57 — . (110) 

m Q J 



C. Local duality violations at m$ ^ 

Here we discuss the impact of new thresholds opening in the spectral density of the vector 
ip^^ip currents, as uiq increases. At = the duality in the (^j^, ^7^) correlation 
function is perfect, since, due to bosonization, the theoretical expression for this correlator 
reduces to exactly one massless state propagation with a known coupling constant. If 7^ 
then higher [ipip] mesonic states appear in the imaginary part of (tp^^ip^lu^) with residues 
proportional to m^. 

To calculate the dependence of the total width Th q on tuq near thresholds, i.e. near 

M Hq = M [ ^ ] + M [ f ap] , (111) 



we need to find the exclusive width T n k(HQ — > [ipip] n + [qq sp }k)- We did it in Sec. |V B| for 
m^p = when only the n = massless <fi state is produced in the [ipip] channel. 

Now, the highly excited states can appear in this channel. They are pseudoscalar ones 
corresponding to even values of n (scalar states are not produced by the conserved current 
ip'y fl ip). Moreover, with nonvanishing Mfif"^ the amplitude of transition to scalar [qq sp \ states 
(odd k) is not proportional to \p\. Thus, near the threshold ( |111| ) the exclusive width of 



decay into the pair of the pseudoscalar [ifjifj] and the scalar [qq sp ] is singular because the factor 
l/|p| in the phase space that explode at thresholds. Therefore, exactly at threshold, T n k 
is infinite. It goes without saying that it cannot coincide with the smooth OPE prediction 
near the threshold. These spikes are clearly visible on the plots of Ref. JTTJ . 

To maximize the exclusive width F n j, for the decay H q — > [if)ip} n + [qq sp \k near the thresh- 
old ( |111[ ) we choose the range of and M^^ close to the corresponding partonic con- 



figuration. In terms of partons we deal with the decay Q — > ip + ip + q. If we fix the square of 
the effective mass of ipip pair 7 then the momentum of the light q is (m^ — / (2t?iq). 

For hadrons close to Mi r and the characteristic hadronic mass squared in qq sp system 

is j ~ /5(mg — M\£)jmQ. Then, it is simple to check that the threshold range is 



M Hq - Mj*® ~ (3 , Mp p] ~ . (112) 

Thus, it is just a few first scalar excitations in qq sp system that are relevant. If Mf® and 
M [ ^ 3p] faUc 
suppressed. 



M^ 9spi fall outside the kinematics indicated above then the exclusive width T n k is additionally 
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Let us remind that the bulk contribution into the total width is provided by the ground 
state in the ipip channel and M^ 3p ^ < ^vtiq3. What is under discussion is a small tail of 

highly exited ipip states which determines the oscillating component. 
The relevant matrix element is 

([#]„fe] fc |(?7^)(W)l^> = ([WWrm ([qq sp }k\qi,Q\ H Q) ■ (U3) 
It is easy to see that the residue ([V>V>] n |'07'V|O) scales as \/W c m^fi j 'M^\ while the transi- 



tion amplitude {[qq S p\k\qiiiQ\HQ) scales as \Jm.Q/3. Assembling all factors together we find 
that near the threshold 

tyi 2 3^ tti 2 3^ 

T nk {H Q - [#]„ + \qq sp ] k ) ~ G 2 N C ^- ~ T Q . (114) 

This equation describes both, the singularity at thresholds and deviation from local duality 
between the thresholds. The spatial momentum \p\ is 



\p\ » pM [ f sp] (M Hq - - M^ p] ) . (115) 

In the middle between the thresholds in the value of \p\ ~ (^/mg) 1 ^ 2 , and 



m 



^3/2 



r n k(H Q -> + [gg sP ]fc) ~ r Q ^ 7/2 . (116) 

This estimate gives the amplitude of the oscillating component. It is applicable also to 
the threshold spikes provided these spikes are averaged over the intervals of M Hq less but 
comparable with the period of oscillations in M Hq dependence (~ tt 2 3 2 /M Hq ). 

Therefore, we conclude that the violation of local duality dies off as 1/toq 2 . This effect 
is by far the largest duality violating contribution. The occurrence of a relatively weak 
suppression is due to (a) the zero resonance width approximation; (b) the singular nature of 
the two-body phase space in two dimensions. Both features have no parallel in actual QCD. 

In summary, we identified two leading effects that are responsible for deviations from 
the parton formula - one is associated with the four-fermion operator in OPE, appearing 
due to m$ ^ 0, the second is the additional duality violating component that was absent at 
rrif = 0. The first, inclusive one, dies off as l/nig, the second (exclusive) at least as l/m 7 q 2 . 
We do not see any room for 1/mg deviations, even oscillating. 



VII. DISCUSSION AND CONCLUSIONS 



The situation we encounter in the 't Hooft model is very instructive. The model is 
readily treatable, which allows one to advance quite far in constructing the OPE series. It 
is superrenormalizable, thus providing an especially clean environment for testing various 
subtle aspects of OPE. The perturbative series for the coefficient functions in the large N c 
limit converges. We find, with satisfaction, that all general statements regarding OPE are 
fully confirmed. 
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The model also clearly exhibits the breaking of local duality by oscillating terms. These 
oscillations are related to the exponential terms in the Euclidian domain and not seen in 
OPE. Due to zero meson widths in the large N c limit they are suppressed only by powers 
of l/rriq which we have determined. 

We note that in the t' Hooft model the local duality of the OPE predictions in the 
inclusive heavy flavor decays (both semileptonic and nonleptonic) holds much better than 
in R(e + e~), the generally recognized classical laboratory for applications of OPE. It is in 
contrast to the opinion often expressed in the literature that OPE is not applicable in the 
inclusive heavy quark widths. Moreover, the numerical computations of Ref. |nj suggest 
that OPE width approaches the (smeared) hadronic ones at a few percent level very soon, 
right after a few first channels are open. 

In actual QCD, already the first excited states are broad enough and inconspicuous, 
leave alone high excitations. When a finite resonance width is introduced, it immediately 
leads to dynamical smearing of the spectral densities, ensuring an exponential suppression 
of the oscillating duality violating terms (see Sec. [VD| and Ref. fL2| ). Thus, in terms of 



actual QCD we are still very far from the solution of this extremely important problem; the 
exercise performed gives us some kind of an upper bound. 

As previously, in actual QCD, we have to rely on models while estimating the expo- 
nential/oscillating terms not seen in OPE. The choice is not large - only two models were 
suggested previously. One of them is an instanton-based model , another is a resonance- 
based model fL2f , close in spirit to estimates we have presented above. The instanton-based 



model is simple and predictive, but it apparently lacks the sophistication inherent to the 
phenomenon in actuality. In particular, it predicts an oscillating component ~ sinmQ, 
rather than ~ sin m 2 Q as would be natural from the resonance point of view. Thus, the 
't Hooft model teaches us that the instanton-based estimates cannot be fully true. On the 
other hand, the resonance-based model, which works satisfactorily in the limit of infinitely 
narrow resonances, does not give a full answer as to how strongly the oscillating component 
is suppressed when the finite resonance widths are switched on. It is clear that further steps 
in developing the existing or engineering new models are needed. 

This work presents the first estimate of the duality violations, from the resonance-related 
considerations based on a 1/N C expansion, in the practically important problem of the 
hadronic r decays. Although not fully conclusive, the results are very encouraging, and call 
for expansion of these ideas in other processes. This is an obvious task for the future. 
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